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The field equations for a time dependent cylindrical cosmic string coupled to gravity are reformu- 
lated in terms of geometrical variables defined on a 2 + 1-dimensional spacetime by using the method 
of Geroch decomposition. Unlike the 4-dimensional spacetime the reduced case is asymptotically 
flat. A numerical method for solving the field equations which involves conformally compactifying 
the space and including null infinity as part of the grid is described. It is shown that the code 
reproduces the results of a number of vacuum solutions with one or two degrees of freedom. In 
the final section the interaction between the cosmic string and a pulse of gravitational radiation is 
■ briefly described. This is fully analyzed in the sequel. 
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I. INTRODUCTION 



Cosmic strings are topological defects that formed during phase transitions in the early universe. They are important 
because they are predicted by grand unified theories and produce density perturbations in the early universe that 
might be important in the formation of galaxies and other large scale structures JjJ. They are also important since 
they are thought to be sources of gravitational radiation due to rapid oscillatory motion Q . In the simplest case of a 
string moving in a fixed background one can take the thin string limit and the dynamics are given by the Nambu-Goto 
action || which is known to admit oscillatory solutions. However in order to fully understand the behavior of cosmic 
, strings one should study the field equations for a cosmic string coupled to Einstein's equations. 

A cosmic string is described by a U(l) gauge vector field coupled to a complex scalar field <E> = \j \[2Se % ^ . The 
Lagrangian for these coupled fields is given by 



7 L M = l-W^SW^S+l-S 2 (W^ + eA ll )(W^cl > + eA^~X(S 2 - V 2 ) 2 -]F^F^, (1) 
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where F^ v — V ' ^A v — V„^4 M , e, A are positive coupling constants and r\ is the vacuum expectation value. 

The Einstein-scalar-gauge field equations for an infinitely long static cosmic string have been investigated by Laguna 
and Garfinkle Q . Some geometrical techniques related to those used in this paper have also been used in the analysis 
of cosmic strings by Peter and Carter j|] and by Carter ||. I n the present paper (and its sequel) we will investigate 
the behavior of a time dependent cylindrical cosmic string coupled to gravity. In particular we investigate the effect 
of a pulse of gravitational radiation on an initially static cosmic string and the corresponding gravitational radiation 
that is emitted as a result of oscillations in the string. Since we are interested in the gravitational radiation produced 
by the string it is desirable to measure this at null infinity where the gravitational flux is unambiguously defined and 
one does not need to impose artificial outgoing radiation conditions at the edge of the numerical grid. However the 
infinite length of the string in the z-direction prevents the spacetime from being asymptotically flat. We therefore 
follow the approach of Clarke et al. [Q] and use a Geroch decomposition with respect to the Killing vector in the 
^-direction to reformulate the problem in 2 + 1 dimensions. Note however that unlike Clarke et al. (?|] we apply the 
Geroch decomposition to the entire problem not just to the exterior characteristic region. The gravitational degrees 
of freedom of the 3 + 1 problem are then encoded in two geometrically defined scalar fields defined on the 2 + 1- 
dimensional spacetime. These are the norm of the Killing vector v and the Geroch potential r for the rotation. We 
show in section III that the energy-momentum tensor of these fields describes the gravitational energy of the original 
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cylindrical problem. An important feature of the reduced 2 + 1 spacetime is that it is asymptotically fiat and this 
allows us to conformally compactify the spacetime and include null infinity as part of the numerical grid. 

In section II we briefly describe the Geroch decomposition and the field equations that one obtains in the cylindrical 
case. We also show how it is possible to rescale t, p and the matter variables to simplify the equations. In order 
to demonstrate the numerical accuracy of the full code it is useful to compare it with either an exact solution or 
else with some other independently produced numerical results. Unfortunately this is not possible for the full code 
where we must satisfy ourselves with the internal consistency of the code as demonstrated by convergence testing for 
example. If one considers only the gravitational part of the code it is possible to specialize to the case of a vacuum 
solution. The simplest of these is the Weber- Wheeler solution ||, but we also consider a rotating vacuum solution 
due to Xanthopoulos || which describes a thin cosmic string and a second rotating vacuum solution due to Piran et 
al. Jio| . In order to compare our results with these exact solutions we must first write them in terms of the Geroch 
potential, this is done in section IV. This description is useful not just for numerical purposes but also in interpreting 
these solutions since the two polarization states for the gravitational radiation have a simple interpretation in terms 
of the Geroch variables v and r. The numerical code used in this case, the convergence analysis and the comparison 
between the numerical results and the exact solution is given in section IV. As far as the matter part of the code is 
concerned there are no exact solutions and one must compare the code with other numerical results. The simplest 
special case involves decoupling the matter variables from the metric variables and considering the equations of motion 
in Minkowski space. There do not seem to be other results available in the dynamic case, but the static solutions 
have been investigated by a number of authors [jy], Jl^] and [jl3). Another solution which has been investigated 
previously is a static string which is coupled to the gravitational field. Finding such solutions is much harder than 
one might suppose due to the asymptotic behavior of the matter variables. As well as the physical solution to the 
equations there is an exponentially diverging non-physical solution [ p"3"[ which must be suppressed. By compactifying 
the radial coordinate we can control the behavior of the solution at infinity in the static case by using a relaxation 
scheme. This allows us to obtain solutions for all values of the radius rather than the fairly restricted range of p 
that had been previously obtained, and also permits us to use the proper boundary conditions at infinity. The static 
string is described in further detail in sections V and VI. Finally we briefly describe the numerics for the dynamic 
string coupled to gravity. Here again the asymptotics of the matter variables make it hard to write a stable code, 
but by using an implicit scheme we are able to produce a code with long term stability and second order convergence 
which agrees with the previous results in the special cases described earlier. From a physical point of view the most 
interesting feature of this code is that it is able to describe the interaction of a gravitational field with two degrees 
of freedom with the full non-linear cosmic string equations. So far we have only investigated the interaction of the 
cosmic string with an incoming Weber- Wheeler type pulse of gravitational radiation with just one degree of freedom. 
We find that the pulse excites the cosmic string and causes the scalar and vector fields to vibrate with a frequency 
which is roughly proportional to their respective masses. This oscillation slowly decays and the string eventually 
returns to its previous static state. This is briefly described in section VII and in detail in the sequel |l4|], which we 
henceforth will refer to as paper II, where comparisons with other results and a full convergence analysis is given. 



II. THE GEROCH DECOMPOSITION 



As we explained above it is not possible for a cylindrically symmetric cosmic string to be asymptotically flat due 
to the infinite extent of the string in the z-direction. By factoring out this direction we can obtain a 3-dimensional 
spacetime that is asymptotically flat. If the Killing vector in the z-direction is hypersurface orthogonal then one can 
simply project onto the surfaces S given by z — const. However we wish to consider cylindrical solutions which also 
have a rotating mode and in this case the Killing vector £ M is not hypersurface orthogonal. Geroch [|5| has shown 
how to factor out the Killing direction in this more general case. The idea is to identify points which lie on the same 
integral curves of the Killing vector field and thus obtain S as a quotient space rather than as a subspace. There is 
then a one-to-one correspondence between tensor fields on S and tensor fields on the 4-dimensional manifold M which 
have vanishing contraction with the Killing vector and also vanishing Lie derivative along £ M . One may therefore use 
the four dimensional metric to define a metric on S according to the equation 

h^ = g^ + {e^r l Uv (2) 
The extra information in the 4-metric is encoded in two new geometric variables; the norm of the Killing vector 

v = -e^ (3) 

(where we have introduced the minus sign to make v positive in the spacelike case) and the twist 
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r M - -e^fVT- (4) 

Geroch then showed how it is possible to rewrite Einstein's equations in terms of the 3-dimensional Ricci curvature 
of (<S, h) and equations involving the 3-dimensional covariant derivatives of v and the twist. If we let define the 
covariant derivative with respect to then one can show that 

D [p r a] = e p ^ v eKC (5) 

where R V T is the 4-dimensional Ricci tensor. It is clear that this vanishes in vacuum so that T a is curl free and may 
be defined in terms of a potential. It is a remarkable fact that this remains true for spacetimes with a cosmic string 
energy-momentum tensor so that even in this case we may write 

T a = D a T, (6) 

where we have introduced the convention of using Latin indices to describe quantities defined on S. 

We may now write Einstein's equations for the 4-dimensional spacetime (M, g) in terms of the Ricci curvature of 
(S, h) and the two scalar fields v and r defined on 5. We obtain 

K ab = ^iy- 2 [(D a r)(D b r) - h ab (D m r)(D m r)} + ^u~ l D a D b v - -v-\D a v){D b v) 

+8nh a ' i h b ' , {T ia) - -g^T), (7) 

D 2 v = ^v- l {D m v){D m v) - u-\D m T){D m T) + 16tt(T^ - ^„T)^f, (8) 

D\=^v- l {D m r){D m v). (9) 

The transformation to the 3-dimensional description is not only mathematically convenient but is physically mean- 
ingful. If the Killing vector is hypersurface orthogonal then r vanishes and the gravitational radiation has only one 
polarisation which may be defined in a simple way in terms of v. If there are both polarisations present then the 
+ mode is given in terms of v while the x mode is given in terms of r [see equation ( |l9| ) below]. One can simplify 
things by making a conformal transformation and using the metric h ab = vh ab in which case we can write the vacuum 
Einstein-Hilbert Lagrangian on M in terms of 3-dimensional variables on S 



Ig= / RV^gd^x (10) 



t {n~\v- 2 [h a \D a T)(D b T) + h a \b aV ){D b v)]}VlSx (11) 
Js 1 



and we see that the 4-dimensional gravitational field is described in three dimensions by the two scalar fields v and r 
conformally coupled to the 3-dimensional spacetime with metric h ab . Since the 3-dimensional spacetime has no Weyl 
curvature it is essentially non-dynamic and we see that v and r encode the two gravitational degrees of freedom in 
the original spacetime. The corresponding 'energy-momentum' tensor for these fields is 

f ab = x - v - 2 \b a rb b T - h ab h cd (b C T){b d T) + b a vb bV - x -h ab k cd (b c v){b d v)\ (12) 

and if there is matter present in four dimensions there are also additional matter terms in three dimensions. As 
shown in equation ( p3| ) this 3-dimensional 'energy-momentum' tensor for v and r gives the correct expression for the 
4-dimensional gravitational energy. 



III. THE FIELD EQUATIONS 



For the case of a cylindrically symmetric vacuum spacetime one can write the metric in Jordan, Ehlers, Kundt and 
Kompaneets (JEKK) form f§, || 

ds 2 = e 2 ^\dt 2 - dp 2 ) - p 2 e^dtf - e 2 ^M0 + dz) 2 . (13) 

However this form of the metric is not compatible with the cosmic string energy momentum tensor so we follow 
Marder ]l8|] by introducing an extra variable p into the metric and writing it in the form 
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ds 2 



^-^(dt 2 -dp 2 ) 



P 2 e 



(14) 



This form of the line element enables us to make easy comparisons with the JEKK vacuum solutions previously 
considered numerically by Dubai et al. JnJ and d'Inverno et al. pQ ]. The field equation for p decouples from those 
for the other metric variables and it has a source term given by Too ~T\\. The physical interpretation of p is briefly 



discussed by Marder 



In the static case one can show that C regularity on the axis implies that p = ln(z/) — 7. 



The metric given by ( |14| ) has zero shift and lapse determined by the condition g u = ~9pp- In this gauge the null 
geodesies are given by the simple conditions u = t — p = const, and v = t + p = const. The remaining coordinate 
freedom is given by the freedom to relabel the radial null surfaces: u — > f(u) and v — > g(v) where / and g are arbitrary 
functions. We may fix this by specifying the initial values of p and its derivative. For example we can choose p to be 
equal to its static value and p^ to vanish, but due to time dependent matter source terms in the evolution equation 
for p [see equation ( p8| ) below] this does not make p constant in time. 

In terms of these variables we find the norm of the Killing vector in the z-direction £^ — <5g to be given by 



and the twist potential is related to 10 by 



D a T = p e 



w p ,w t ,0,0). 



Finally the conformal 3-metric h a b is given by 

da 2 = e 2 ^(dt 2 - dp 2 ) - p 2 e 2 »d<f> 2 . 
It is also of interest to calculate T a f, in terms of these variables. We find 



T a bt a t 



a ib 



1 



-2(7+M) 



where t a is a unit timelike vector proportional to Jj. Note that the quantities 



A = 



B = 



c = 



D=(^ 
v 



(15) 



(16) 



(17) 



(18) 



(19) 



are exactly the same as the quantities A, B, C and D which are given (by more complicated expressions) in terms 
of ip an d u) in equations (4a)-(4d) of Piran et al. |^IJ and describe the two polarizations of the cylindrical gravita- 
tional field. Furthermore if we consider the special case of vacuum solutions and integrate T a f 1 t a t b over the region 
V = {0 < p < po, t = to} with respect to the volume form dV on t — to we find 



E{t ,p ) 



T ab t a t b dV 



7T 

4 

2tt 



v 

Pa 



"f,P e 



-~<dp 



er+(^r+(i7 



pdp 



(20) 

(21) 

(22) 
(23) 



where we have used the vacuum field equations for 7^ [sec equation (|3lj) below]. Note that this is the same as the 
energy obtained by Ashtekar et al. |^2| but does not require the Killing vector to be hypersurface orthogonal. It differs 
from the C-energy in general but agrees with it in the linearized case. 

As far as the matter variables are concerned we make the obvious generalization of the form used by Garfinkle |26] ] 
and write 



A„ = V(i,p)-l]V, 



(24) 
(25) 



We may now write the field equations for the complete system. Since we are working in three dimensions we have 
three independent evolution equations. After some algebra these may be written as 



4 



-STr^Xu-'e^+^iS 2 - rf) 2 + e- 2 p- 2 ve-^{P\ - P 2 p )], 



Up, + p"V,p - fJ- 2 t + M 2 P = -8tt[2Ai 



0, 



-l 2(7+^/02 



{S 2 ~r ] 2 f+p- 2 e^S 2 P 2 }, 



where □ represents the flat spacetime d'Alembertian which in cylindrical coordinates is given by 

cP_ <P_ _ x d_ 
dt 2 + d P 2+p dp' 



□ 



(26) 
(27) 
(28) 



(29) 



There are also three constraint equations, one of which vanishes identically due to the rotational symmetry The 
remaining equations give 



7,t 
7, P 



1 + PM„ 



P 2 A ~ (1 + PM,p) 2 



(J>,t P - P>,t{l,p + M,p) + \v 2 (~r,tT,p + v,tv^p) + 8iT(S,tS,p + e 2 p 2 ve 2fl P t P p ) 
((1 + pp, p ) { - 4 7 r[2,- 1 e 2 ^^)A(5 2 - V 2 ) 2 + (S 2 t + S 2 p ) 



(30) 



+e- 2 p- 2 ve- 2 ^P 2 + P 2 p ) + p- 2 e 2 ^S 2 P 2 ] ^(r 2 + r 2 p + v 2 + v 2 p ) p~^, p - p, pp 



2^ 2 ( T ,tT. p + v tt v, p ) + p V,t + ^,tp + 8n(S tt S,p 



e p ve 



-2p 



(31) 



Finally there are the equations for the matter variables S and P which may be derived from the Euler-Lagrange 
equations for Lm or alternatively from the contracted Bianchi identities. 



US - S,tli, t + S <p p >p = SiAX^e^+^iS 2 - T) 2 ) + p- 2 e 2 ->P% 

UP - P, t {v-\t ~ /x,t) + PA^'^p - »,p - 2 P _1 ) = e 2 is- 1 e 2 ^+^PS 2 



(32) 
(33) 



We also need to supplement these equation by boundary conditions on the axis. For the 4-dimensional metric 
variables the simplest condition is to require the metric to be C 2 on the axis so that we have a well defined curvature 
tensor. This gives the conditions 



In terms of v and r this gives 



tl>{t,p) =ai(i) + <3(p 2 ), 
u;(t,p) = 0(p 2 ), 
p(t,p)=a 2 (t) + 0(p 2 ), 
1 {t,p)=0{p). 



v(t,p) = a 3 (t) + 0(p 2 ), 
r(t,p) = 0(p 2 ), 



(34) 
(35) 
(36) 
(37) 



(38) 
(39) 



where we have chosen the additive constant in the definition of the potential r so that it vanishes on the axis. In 
certain situations C 2 regularity is too strong and one must impose the weaker condition of elementary flatness p3| . 
However even this is too strong for the Xanthopoulos solution which has a conical singularity on the axis. 
The boundary conditions for S and P on the axis are |2q] 



S(t,p) = 0(p), 
P{t,p) = l + 0(p 2 ). 



(40) 
(41) 



In order to consider the behavior of the solution at null infinity we first transform from t to a null time coordinate 
u = t — p. We also wish to compactify the region and following Clarke et al. (t]] we make the transformation 



1 



(42) 
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In terms of these variables the equations become 

= 8M-2X^e^\S 2 - V 2 ) 2 + le- 2 y 7 ve-^{\y*Pl + 2P„P y )], (43) 

□r - y 3 v~ x (r u v, y + T tV v, u + ^y\ v i>, y ) + ^y 3 {fJ,, y T yU + fjt, u r, y + ^y 3 ^, y T y ) = 0, (44) 

= -87r[2Ai/- 1 e 2 (T + '*)(5 2 - ?? 2 ) 2 + y 4 e 2 ^ 2 P 2 ], (45) 



-8tt{(1 - \yii >y )[S% + iy 3 5 u S, y + e" 2 */ 4 wT 2 "(P 2 + -y^P,)] 

-\v\u{Sl + ^V^P 2 )})/^ 2 - ( M ,„ + \^ y )] 2 - /xfj, (46) 
7,«= [- J(32/ 2 M,, + yV, w ) + \v 3 V 2 y - \y*v- 2 {rl + v 2 y ) - 2ny 3 (S 2 y + e" 2 */ 4 wT 2 "P 2 )} 

/(y 2 ~\y 3 ^ y ), (47) 

OS + \y 3 S, u ^ y + \y 3 S. y {^ u + l -y 3 ^ y ) = S^Xv' 1 e 2 ^{S 2 - V 2 ) + y^P 2 }., (48) 

OP + P,u[\v 3 {^\v ' + V] + ^PyW'Hvu + \y 3 ^y) - + \y 3 ^,y) + ^} 
= e 2 v~ 1 e 2 ^PS 2 , (49) 
where in these coordinates □ is given by 

D=y i{ 4y i&y +y4 $ +y3 Yy~ 4 i)- (50) 

It is worth remarking that one can have solutions to these equations which are regular at y — and which represent 
4-dimensional metrics in which u> diverges. Thus the notion of the 3-dimensional spacetime being asymptotically flat 
is weaker than might first be supposed. The asymptotic behavior of S and P at null infinity is given by 

S(u,y)=r) + 0(y), (51) 
P(u,y)=0(y). (52) 

These are discussed in more detail in section V. 

The field equations are solved numerically in two ways. Firstly an explicit second order Cauchy-Characteristic 
Matching (CCM) scheme similar to that employed by Dubai et al. |l9| and d'Inverno et al. p(| is used, but using 
a Geroch decomposition in the whole spacetime (not just the characteristic portion) which allows one to use the 
geometrically defined variables in both the interior and exterior regions. This scheme works very well when compared 
to the exact vacuum solutions but is less satisfactory when the matter terms are included (see below). An alternative 
scheme with similar accuracy but with long term stability is a fully characteristic second order implicit scheme. This 
has the advantage that the scheme naturally controls the growth of the derivatives at infinity and hence automatically 
selects the physical rather than the non-physical solutions. The details of this scheme are described in paper II. 

IV. EXACT VACUUM SOLUTIONS 

In this section we describe the exact vacuum solutions which are used to test the codes. The solutions we will 
consider are the Weber- Wheeler gravitational wave || which just has the + polarization mode and two solutions due 
to Xanthopoulos || and Piran et al. (T^j which have both the + and x polarization mode. 



G 



The first exact solution we consider is the Weber- Wheeler gravitational wave originally investigated by Einstein and 
Rosen [^4| . It consists of a cylindrically symmetric vacuum wave with one radiational degree of freedom corresponding 
to the + polarization mode M] . It describes a gravitational pulse originating from past null infinity and moving toward 
the z-axis. After imploding on the axis, it emanates to future null infinity. 

This solution has no rotation so that u> and hence r vanish and the solution may be described in terms of tp which 
satisfies the wave equation. A solution to the wave equation in cylindrical coordinates may be given in terms of Bessel 
functions and by superposing such solutions we may write 



t/j{t,p) = 2b I e- an J (np)cos(nt)dn, 



(53) 



where a > 0. For convenience we let 



X = a 2 + p 2 — t 2 



and one may show that (p3) may be written in the alternative form 



u(t, p) = exp 



2b\ 



1 2{X + ^ X 2 + 4a 2 t 2 ) 
X 2 + 4a 2 t 2 



(54) 



(55) 



The corresponding value of 7 is obtained by integrating 7 iP = p(ip 2 t + tp p ) and using jit, 0) — and is found to be 

X 2 - Aa 2 t 2 



7(t,p) 



2a 2 



1 - 2a 2 p 2 



a 2 + t 2 - p 2 



(X 2 + 4a 2 t 2 ) 2 VA 2 + 4a 2 t 2 



(56) 



The next solution we consider is one due to Xanthopoulos [Q] which has a conical singularity on the z-axis and 
therefore describes a rotating vacuum solution with a cosmic string type singularity. Xanthopoulos derived the 
spacetime by finding a solution to the Ernst equation in prolate spheroidal coordinates. To compare this with our 
numerical result we must transform to cylindrical coordinates and also find the Geroch potential. It is convenient to 
first define the following quantities 



Q = p 2 -t 2 



1, 



X = + 



Y = -[{2a 2 + l)X + Q] + l- a^2(X-Q), 



Z 



-[(2a 2 + l)A + Q]-l, 



where < \a\ < 00. The solution derived by Xanthopoulos then becomes 

,, s 1 , Z 



u(t, p) 

l(t,p) 



y/a 2 TT{X + Q-2){Z-Y) 



2aZ 



— In — 7; — , 

2 a 2 X : 



(57) 
(58) 

(59) 
(60) 



(61) 
(62) 
(63) 



where we have imposed j(t, 0) = 0. A straightforward but rather tedious calculation shows that this satisfies Einstein's 
field equations. (This and a number of other calculations in the paper were checked using the algebraic computing 



package GRTensor II 25 



The norm of the Killing vector in the z-direction is given by 

v(t,p) 

The Geroch potential is easily obtained from the Ernst potential and is found to be 



Z 
Y' 



(64) 
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r(t,p) 



v/2(a 2 + 1)VXTQ 



Y 



Note that v and r satisfy (E6[)— rt2 



m vacuum, i.e. 

,-l/_2 



r— 1 —11 1 Z Z 



= 0, 



□t + 2i/ L (T, t u, t - T tP U tP ) = 0. 



(65) 



(66) 
(67) 



Expressions for all these quantities may be obtained in the exterior characteristic region by transforming to the 
(it, y) variables. Although ip tends to zero as one approaches null infinity, to diverges so that the 4-dimensional metric 
is not asymptotically flat even along null geodesies lying in the planes z = const. However by contrast the Geroch 
potential r vanishes as one approaches null infinity in the 3-dimensional spacetime. This is an example of the fact 
that the Geroch potential can be well behaved even if oj in the JEKK form of the 4-metric diverges as one goes 
outward in a null direction. We also give an expression for the gravitational flux at infinity which is given by E tU 



where E{u, y) = 27r[l 



>] 



lim E u = 



4tt 



(1 + u 2 )[(l + 2a 2 )\ZT~H? - u] 



< 0. 



(68) 



Thus the string is losing energy through gravitational radiation. 
To plot the solution for < p < oo we introduce the radial variable 



VVp 



for < p < 1 
for p > 1 , 



(69) 



thus < w < 3 where the infinite value of p is mapped to w = 3. This choice is slightly different from that of Dubai 
et al. |l9| and avoids discontinuities in the radial derivatives at the interface due to the square root in the definition 
of y. Plots of v, t and 7 as given by (|63|)-(65) are shown in Fig. 1. The error in the numerical results as computed 
using the CCM code are shown in Fig. 4. 





FIG. 1. The exact Xanthopoulos solution for —70 < t < 30, < w < 3 and a = 0.5. Plots are from left to right: u(t,w), 
r(t, w) and y(t, w). One clearly sees the incoming pulses in v and r. As the pulse hits the axis the string loses energy through 
gravitational radiation as seen in 7. 

The final exact solution we consider is one due to Piran et al. [[To) which also has two degrees of freedom representing 
the two polarization states. As in the case of Xanthopoulos' solution, it represents two incoming pulses that implode 
on the axis and then move away from it. Piran et al. obtained their solution by starting with the Kerr metric in 
Boyer-Linquist form, transforming to cylindrical polar coordinates and then swapping the t and z coordinates (and 
introducing some factors of i to maintain a real Lorentzian metric). See JlCj ] for details. The resulting metric may 
be written in JEKK form. The solution is rather complicated but may be simplified by introducing the following 
additional quantities 



R = b- 1 W& + (t- p) 2 -t + p], 

s = b- 1 Wb 2 + (t + P ) 2 + t + p], 



T = 1 + RS + 2a~V(a 2 = 1)^ 



(70) 
(71) 
(72) 



and 
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X={l + R 2 ){l + S 2 ), 
Y = a 2 T 2 + (R-S) 2 , 
Z = a 2 (l~RS) 2 + {R + S) 2 , 

where 1 < a < oo and < b < oo. The metric coefficients are then given by 



,/ \ 1 , z 
V\t,p) = 2 ln y' 



w(t,p) = by/ a 2 - 1 



7(*^)=2 ln x- 



2 1 



1^ (R + S) 2 T 

Vrsz 



(73) 
(74) 
(75) 



(76) 
(77) 
(78) 



Notice that Minkowski space is obtained in the limit that a — ► 1, and that we can also consider the case b — ► in 
which case the rotation vanishes. This is not true for the Xanthopoulos solution which is only real for a sufficiently 
large rotation. 




The solution is regular on the axis, but like the Xanthopoulos solution to diverges as one approaches null infinity. 
Again the answer is to transform to the v, r variables which are regular both on the axis and at null infinity. Finding v 
is straightforward; however solving the differential equations for the Geroch potential t for such a complicated metric 
is extremely difficult, but r may be found by first finding the Geroch potential for the timelike Killing vector of the 
Kerr solution and then making the appropriate transformations. Note that the same process transforms the Killing 
vector into one along the z-axis. One then finds 

"(t,p) = y, (79) 



r(t, P )=- , . (80) 

[2v/(a 2 - 1)RS + o(l + RS)} 2 + (R - S) 2 

The corresponding results in the characteristic region are easily found by transforming to (it, y) coordinates. Plots of 
v, t and 7 as given by ( |7^ ) (jsfj) are shown in Fig. 2. The error in the numerical results as computed using the CCM 
code are shown in Fig. 5. 

We now briefly describe the accuracy and the convergence analysis for the explicit CCM version of our code. The 
results for the the implicit version are similar and are given in paper II. We define the pointwise error at the ith time 
slice and jth grid point for some function / by 

£j(f) = f{U, Wj) eX act - f(ti, Wj) computed- (81) 

The pointwise error for the vacuum solutions is shown in Fig. 3-5 for 600 grid points and 10,000 time steps corre- 
sponding to < t < 15. The code is stable and accurate for at least 20,000 time steps with a Courant factor of 0.45. 
Beyond this point the metric functions have almost decayed to zero and the dynamical behaviour is very slow. 
In order to analyze the convergence of the code we define the spacetime i^-norm for some function / as 



JVi-JV 2 



nm=\ rTz n , (82) 
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FIG. 3. Pointwise error for the Weber-Wheeler solution. From left to right: tp(t,w) ■ 10 6 and j(t,w) ■ 10 7 . 




FIG. 4. Pointwise error for the Xanthopoulos solution. From left to right: v(t,w) ■ 10 5 and r(t,w) ■ 10 5 . 



where N x is the number of time slices and N 2 is the number of grid points on each slice. We also define the relative 
norm 



fr = 



^2i : j[f j exact! 



(83) 



Convergence testing of the code is done by doubling the grid size keeping the Courant factor constant. We therefore 
also need to double the number of time steps. We measure the convergence through 



(84) 



For a convergent second order code one should have f c — 4, and we can clearly see from Tables I— III that the code 
has achieved second order convergence in time and space. 



Grid pts. 


v T 


7r 


Time steps 


v c 


7c 


300 


6.29 


10" 


-Y 


9.38-10" 


-6 


5,000 






600 


1.11 


10" 


-7 


1.66- 10" 


-6 


10,000 


4.01 


4.01 


1,200 


1.96 


10" 


-8 


2.94- 10" 


-7 


20,000 


4.00 


4.00 


2,400 


3.47 


10" 


-9 


5.19-10" 


-8 


40,000 


4.00 


4.00 


4,800 


6.14- 


io- 


10 


9.18-10" 


-9 


80,000 


4.00 


4.00 



TABLE I. Convergence test 



the Weber- Wheeler solution. 



Grid pts. 


V 






7r 


V c 




7c 


300 


1.70- 


10" b 


5.71 • 10 


-b 


1.66-10 


-b 








600 


2.99- 


io- 7 


1.01 • 10 


-6 


2.93-10 


-7 


4.01 


4.00 


4.01 


1,200 


5.28- 


io- 8 


1.78-10 


-7 


5.18-10 


-8 


4.01 


4.00 


4.01 


2,400 


9.32- 


io- 9 


3.15-10 


-8 


9.16-10 


-9 


4.00 


4.00 


4.00 


4,800 


1.65 • 


io- 9 


5.57-10 


-9 


1.62- 10 


-9 


4.00 


4.00 


4.00 



TABLE II. Convergence test: the Xanthopoulos solution. 
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FIG. 5. Pointwise error for the Piran et al. solution. From left to right: v(t, w) ■ 10 and r(t,w) ■ 10 5 . 



Grid pts. 


V 


r 




T 






7 


r 


v c 




7c 


300 


1.07- 


io- 


b 


3.94- 


10" 


(> 


2.76- 


10~ b 








600 


1.89- 


10" 


6 


6.96- 


io- 


-7 


4.87- 


lO" 7 


4.00 


4.00 


4.01 


1,200 


3.35- 


10" 


7 


1.23 ■ 


io- 


-7 


8.61 • 


io~ 8 


4.00 


4.00 


4.00 


2,400 


5.92- 


10" 


8 


2.17- 


io- 


-8 


1.52- 


io- 8 


4.00 


4.00 


4.00 


4,800 


1.05- 


10" 


8 


3.84- 


io- 


-9 


2.69- 


io- 9 


4.00 


4.00 


4.00 



TABLE III. Convergence test: the Piran et al. solution. 



V. THE COSMIC STRING IN MINKOWSKI SPACETIME 



In this section we examine the field equations for the cosmic string. The simplest case is to look at the equations 
of motion on a fixed Minkowski background. If we do this then the Euler-Lagrange equations for ([!]) give 



US = 5[4A(S* 2 - rf) - p- 2 P\ 
UP - 2p- 1 P p = e 2 S 2 P, 



(85) 
(86) 



where □ is the d'Alembertian in cylindrical polar coordinates given by (p9|). 

It will turn out that the solutions of this simpler set of equations are qualitatively similar to those of the full system 
of equations for a dynamic cosmic string coupled to Einstein's equations. However the full system enables one to 
perturb a static string with a pulse of gravitational radiation and in turn look at the effect of the string's oscillations 
on the gravitational waves. 

A special case of ([35]), (^6|) is when one looks for a static solution. This has been looked at before by a number of 
authors, for example Garfinklc ^6|, Laguna et al. pJ| , Laguna-Castillo et al. [Q and Dyer et al. [^3|. For a static 
string equations (p5|) and (|86|) reduce to 



dp 
d 

^ dp 



pi 



(87) 



dp 



This pair of coupled second order equations requires four boundary conditions, 
energy solution these are 



5(0) = 0, 
P(0) = 1, 



lim S(p) = f], 
lim P(p) = 0. 



For the physically relevant finite 

(89) 
(90) 



It is not possible to obtain an exact solution to these equations but one can investigate the asymptotic behavior 
for large p. The solution satisfying the above boundary conditions has asymptotic behaviour given by |27|] 



S(p)~Ti-Ko(V8\np)~Ti 
P(p) ~ pKxterjp) r 



p 1 e 



/32Ar/ 



er)p 



p-^e- 



/SXrip 



(91) 
(92) 
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However as well as these physical solutions, the equations admit non-physical solutions which have exponentially 
divergent behavior as p — > oo. It is the existence of these non-physical solutions which make the problem rather 
delicate from a numerical point of view and makes a method such as shooting hard to apply. 

Before proceeding further we follow Garfinkle [^6| by introducing rescaled variables and constants which simplify 
the algebra (and are also important when considering the thin string limit). Provided we rescale the time coordinate 
this also simplifies the fully coupled system. Let us introduce 

X = ~, (93) 
V 

r = V\r)p, (94) 
t = y/Xrft, (95) 
e 2 

a = — . (96) 
A 

Thus a represents the relative strength of the coupling between the scalar and vector field given by e compared to 
the self-coupling of the scalar field given by A. Critical coupling, when the masses of the scalar and vector fields are 
equal, is given by a — 8. In the theory of superconductivity, a = 8 corresponds to the interface between type I and 
type II behaviour J27J . 

With the above rescaling equations ( p7|) and ( pq ) become 

's(rf)-^(^-D+n vm 

T L (r-'i£) =aX 2 P. (98) 
dr \ dr J 

Note the rescaled version of (p6|)-(p3|) may be found simply by making the replacements 

S -> X, (99) 

p -> r, (100) 

t -> t, (101) 

e 2 -> a, (102) 

A->1. (103) 

The boundary conditions of the cosmic string are given by 

X(0) = 0, lim X(r) = 1, (104) 

r — »-oo 

P(0) = 1, lim P{r) = 0. (105) 

r — 'OC 

Because of the boundary conditions at infinity it is desirable to introduce a new coordinate which brings in infinity 
to a finite coordinate value. For this purpose we again introduce an inner region (r < 1) and an outer region (r > 1) 
where we use the radial coordinate y given by 

V = 4r- (106) 

In order to solve equations (^), ( |98| ) numerically we introduce a spatial grid consisting of n\ points in the inner 
region and n-i points in the outer region. The points n, . . . , r ni cover the range < r < 1, and y ni +i, ■ ■ ■ ,Vni+n 2 
cover the range 1 > y > 0. Thus, r = 1 = y is represented by two points. These two points form the interface of the 
code where r-derivatives are transformed into y-derivatives. The static equations in the outer region are 



d ( dX 

V dy \ dy 



(X 2 - 1) „ 

4^ + P 2 

V 



(107) 



y±(y* d £)=AaX*P. (108) 

In order to apply boundary conditions at both r — and y = the equations were solved numerically using a 
relaxation scheme (as described in pM] for example). For this purpose we wrote the equations as a first order system 
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in both regions (see paper II) and used second order centered finite differencing. Solutions were computed in this way 
for different resolutions N — m + n-2 to check the convergence of the code. Since there is no exact solution available 
the convergence was checked by calculating the ^2-norm with respect to a high resolution result for N = 2400 (1200 
points in each region) 



N 



(■2 [A/ 



Ah 



2400N 2 



(109) 



where / stands for either X or P. For a second order convergent code one would expect the £2-norm to decrease by 
a factor of 4 if the number of grid points is doubled. We find that the code shows clear second order convergence. 
Since the convergence in this case is very similar to that of a static string coupled to gravity considered in the next 
section, we only give the results for the latter case in Table IV. 




FIG. 6. Cosmic string in Minkowski spacetime. X and P for a: 0.125, 1, 8, 64. As a increases, both X and P become more 
concentrated towards the origin. 

In Fig. 6. we plot X and P for a = 0.125, 1, 8, 64. The results show that for fixed values of the self-coupling A and 
vacuum expectation value r\ of the scalar field, both the vector and scalar fields become more concentrated towards 
the origin as the coupling between the fields e increases. One also finds that for a fixed ratio of the coupling constants 
a the fields become more concentrated towards the origin as either the self-coupling A or the vacuum expectation 
value T) of the scalar field increases. 



VI. THE STATIC COSMIC STRING COUPLED TO GRAVITY 



The next class of solutions we wish to consider are static solutions of the fully coupled equations. In the case of no 
t dependence equations (E6h-(B3h reduce to 



2 2 
V r ~ T r 2 
(rv^ r ). r ~ —V.rP'.r H ' : — h Snrj 



V 



2e 2( 7+Al)(x 2 _ 1)2 



~(»T>),r= V ( 2 ~7 > 



V 

2 Q„„2 



-o{r H,r),r= -M r ~ 87n 7 



,27 



X 2p2 



=>2(7+p) 



1 + rfi ir I Av 



— (rX r ) r = —X r jl r + X 

r ' ' 



( T ,r + v2 r) + 47r? 7 2 



-(x 2 -iy 



p2 x 2 P 2 pZh+v) 
X 2 , + ve- 2 ^ - e 2 ^^4- - 2- (X 2 - l) 2 



A- (X 2 -l) + e 2 ^ 



(110) 

(111) 
(112) 

U,r, (H3) 

(114) 
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P, 



(/V 

V v 



=2(7+M) 



-PX 1 



(115) 



Notice that ( |lll[ ) with the corresponding boundary conditions is satisfied by the trivial solution r = 0. One also has 
from the field equations 



(f*7,r),r = — H>/V + /V + 87T?7' 



2 I e 27 



(116) 



This equation is a direct consequence of the other equations and will not be used in the calculations but is instead 
used as a check for the code. The equations in the outer region as well as the resulting first order system, the interface 
equations and boundary conditions are given in paper II. We use the same grid and numerical method as in the 
Minkowskian case. In order to check the code for convergence we again compute the ^2-norrn with respect to a high 
resolution calculation. The results are shown in Table IV for a = 1 and the large value r\ — 0.1 and clearly indicate 
second order convergence. Small deviations from a convergence factor of 4 are expected since we compare against a 
high resolution reference solution rather than an exact solution. The same result has been obtained for other choices 
of a and ?/. 





V 


i 1 


7 


X 


P 




1.28- 10~ Y 


2.51 • 10-" 


2.39- 10~ b 


4.16- 10- y 


5.95- 10~ Y 




3.56 


3.59 


3.58 


3.37 


4.04 




3.76 


3.79 


3.78 


3.60 


4.19 




4.58 


4.61 


4.60 


4.44 


4.98 



TABLE IV. Convergence test: static cosmic string coupled to gravity. 




FIG. 7. Cosmic string coupled to gravity, a = 1 and r\ — 10 3 . 
multiplied by 10 5 . On the right X and P are plotted. 



On the left the metric variables: j, u — 1 and /i are plotted 



The metric and matter variables for a — 1 and rj — 10~ 3 are shown in Fig. 7. We find that the behaviour of X and 
P is very close to that obtained for a static cosmic string in Minkowski spacetime shown in Fig. 6. For physically 
realistic values of r\ the metric variables at infinity are close to their Minkowskian values although the non-zero value 
°f 7o = linv_ 7(r) indicates that the spacetime is asymptotically conical, that is Minkowski spacetime minus a 
wedge with deficit angle Atf> = 2n(\ — e -70 ). Thus a string with a = 1 and -q = 10 -3 has an angular deficit of about 
2 x 10~ 5 which corresponds to a grand unified symmetry breaking scale of about 10 16 GeV j27[]. For larger values of 
rj, however, the deviation from the Minkowskian case becomes substantial. Close to critical coupling the deficit angle 
exceeds 27r for values of rj greater than about 0.28 which explains why the code converges well for 77 < 0.28. 



VII. THE DYNAMIC COSMIC STRING COUPLED TO GRAVITY 



In this section we discuss the interaction between the cosmic string and the gravitational field. Here we will simply 
outline the numerical scheme, the full details are given in paper II. The field equations in the (i, p) coordinates are 
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given by equations (|26|)-(|33|) while those in (u, y) coordinates are given by equations (^3|)-(|49|). In both cases there 
are two additional Einstein equations which are consequences of the others and are used to check the code. In fact two 
numerical schemes were employed. The first was an explicit CCM scheme similar to that employed by Dubai et al. C3] 
and d'Inverno et al. (20). However the use of the geometrical variables v and r in both the interior Cauchy region and 
the exterior characteristic region significantly improves the interface and results in a genuinely second order scheme 
with good accuracy even with both polarizations present. For the vacuum equations the code also exhibits long term 
stability. However when the matter variables are included this code performs less satisfactorily. This is because of 
the existence of exponentially growing non-physical solutions. It is possible to control these diverging solutions by 
multiplying the it-derivatives of the matter variables by a smooth 'bump function' which vanishes at y — but is 
equal to 1 for y > c (where c is a parameter). This produces satisfactory results but the bump function introduces 
some noise into the scheme which eventually gives rise to instabilities. 

A much better solution is to control the asymptotic behavior at infinity by using an implicit scheme. The main 
problem with the system of differential equations is the irregularity of the equations at both the origin and null 
infinity. Therefore just as in the static case considered in the previous section the scheme employed divides the 
spacetime into two regions, an inner region r < 1 where coordinates (u, r) are used and an exterior region r > 1 where 
the (u, y) coordinates are used. The equations in both regions are written as a first order system connected by an 
interface and the evolution is accomplished using a code based on the implicit Crank-Nicholson scheme. This implicit 
scheme provided a simple way of implementing the boundary conditions and thus circumventing all problems with the 
irregularities. The outer boundary conditions as well as the equations in the outer region, the first order system used 
for the numerics and the interface are given in paper II. The implicit code showed very good agreement with both the 
exact (vacuum) and previously obtained (static) numerical solutions. It also showed clear second order convergence 
and very long term stability. For convenience characteristic coordinates were also used in the inner region but we do 
not think that their use was responsible for the good features of this code and believe that an implicit CCM scheme 
would have produced similar accuracy, convergence and long term stability. 

The code is able to consider the interaction of the cosmic string with a gravitational field with two degrees of 
freedom, however here we simply describe the interaction with a Weber-Wheeler wave which has just one degree of 
freedom. We consider a pulse which comes in from past null infinity and interacts with a cosmic string in its static 
equilibrium configuration. This interaction causes the string to oscillate which in turn affects the gravitational field 
as measured by v and t. The oscillations in both X and P decay as one approaches null infinity (i.e. as r — > oo for 
fixed u) and also for fixed r as u — > oo. After the oscillation has died away the string variables X and P return to 
their static values. Note however that this decay is rather slow and being able to show this effect depends upon the 
long term stability of the code. Although oscillations are observed in both X and P the character and frequency of 
these oscillations is rather different. 




0.5 1 1.5 2 2.5 3 -5 5 10 15 



FIG. 8. On the left: P(u, w) for a = 1, rj = 10 3 and u — 8.5. On the right: The corresponding contour plot of P(u,r) for 
< r < 200 and —8 < u < 16. Note the use of r for values greater than one in this plot. 

In Fig. 8. we plot P for a = 1 and rj = 10~ 3 at a time u = 8.5 (left panel). The oscillations out to large radii can 
be clearly seen . The contour plot on the right shows the ringing behaviour of the string and the slow decay of the 
oscillations in P. In contrast the oscillations in X displayed in Fig. 9. are restricted to small radii and decay on a 
shorter timescale. 

An investigation of the frequencies fx and fp of the oscillations of the scalar field X and the vector field P (in the 
rescaled unphysical variables) indicates that they are relatively insensitive to the value of rj and the Weber- Wheeler 
pulse which excites the string. They are also largely independent of the radius at which they are measured. However, 
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although fx is also independent of a, we find that fp is proportional to \fot. When we convert to the physical fields 
S and P and use the physical coordinates (t, p) we find that the frequencies in natural units are given by 

fs ~ V\V (117) 
fp ~ ev- (118) 

If we now use the fact that the masses of the scalar and vector fields are given by m| ~ A77 2 and mp ~ e 2 T] 2 E^] this 
gives 

fs ~ m s (119) 

fp ~ m P . (120) 

In fact these frequencies are also obtained by considering the simpler model of a dynamic cosmic string in Minkowski 
space with equations of motion (76)-(77) provided one gives S and P similar initial conditions to that produced by 
the interaction with the pulse of gravitational radiation. Thus the main role of the gravitational field as far as the 
string is concerned is to provide a mechanism for perturbing the string. This is discussed more fully in paper II. 



VIII. CONCLUSION 



In this paper we have shown how the method of Geroch decomposition may be used to recast the field equations for 
a time dependent cylindrical cosmic string in four dimensions in terms of fields on a reduced 3-dimensional spacetime. 
This has the advantage that it has a well defined notion of null infinity. It is therefore possible to conformally 
compactify the 3-dimensional spacetime and avoid the need for artificial outgoing radiation conditions. An additional 
feature of this approach is that it naturally introduces two geometrically defined variables v and r which encode the 
two gravitational degrees of freedom in the original spacetime. 

We have described how the field equations for a cosmic string coupled to gravity may be solved numerically by 
dividing the 3-dimensional spacetime into two regions; an interior region r < 1, and an exterior region r > 1 in which 
the coordinate y is used. The use of the geometric variables v and r greatly simplifies the transmission of information 
at the interface r = 1. Although an explicit CCM code worked very effectively in the vacuum case, the asymptotic 
behavior of the matter fields S and P made it less effective when the string is coupled to the gravitational field. An 
alternative implicit fully characteristic scheme however exhibited good accuracy and long term stability. 

In this paper we have demonstrated the effectiveness of the CCM codes in reproducing the results of the Weber- 
Wheeler solution and two vacuum spacetimes with two degrees of freedom due to Xanthopoulos and Piran et al. This 
involved calculating the Geroch potential for these solutions and using it to compare with the numerically computed 
values. We have also described a code for a static cosmic string which uses a relaxation scheme and provides initial 
data for the dynamic code. 

In the final section of the paper we briefly described the interaction of the cosmic string with a Weber- Wheeler 
gravitational wave. The pulse of gravitational radiation excites the string and causes the fields S and P to oscillate 
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with frequencies proportional to their respective masses. The full details of the code, the convergence testing and the 
interaction between the string and the gravitational field are described in paper II. 
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